Dynamic Hubbard Model: Effect of Finite Boson Frequency 
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Dynamic Hubbard models describe coupling of a boson degree of freedom to the on-site electronic 
double occupancy. In the limit of infinite boson frequency this coupling gives rise to a correlated 
hopping term in the effective Hamiltonian and to superconductivity when the Fermi level is near 
the top of the band. Here we study the effect of finite boson frequency through a generalized Lang- 
Firsov transformation and a high frequency expansion. It is found that finite frequency enhances 
the tendency to superconductivity in this model. 
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INTRODUCTION 

An attempt to understand electron-electron interac- 
tions in solids over the last few decades has focused on 
simplified models which single out, for example, on-site 
or nearest neighbour Coulombic interactions. While it 
is usually understood that the parameters that are re- 
quired in these models are not 'bare' parameters, it is 
often tacitly assumed that the omitted physics can be 
recovered by simply modifying these parameters to ef- 
fectively incorporate the missing physics. However, the 
notion that this procedure may in fact leave out much 
of the interesting physics has been emphasized by one of 
us recently and a new class of model Hamiltonians, 
'dynamic Hubbard models' , was introduced to remedy 
this difficulty P, 0,0,0] . These Hamiltonians describe an 
essential difference between the empty and the doubly oc- 
cupied Wannier orbital, that the former is representable 
by a single Slater determinant and the latter one is not. 
This difference is not described by the conventional mod- 
els and leads to electron-hole asymmetry and to an un- 
coventional mechanism of superconductivity 

In this paper we focus attention on one particular dy- 
namic Hubbard model, which introduces a modulation 
of the Hubbard U by coupling the electron double occu- 
pancy to a fictitious local boson displacement. The site 
Hamiltonian that describes this model is given by Q: 
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where the boson is characterized by an Einstein oscillator 
with mass M and spring constant K. The displacement 
of the local oscillator qi is coupled to the electron dou- 
ble occupancy with coupling constant a, and as a conse- 
quence the on-site repulsion U +aqt becomes a dynamical 
variable . The Hamiltonian in electron representation is 
given by 



H = ui 
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where cj CT is an electron creation operator for site i and 
spin (j, riiu is the corresponding number operator, and 
we have introduced the boson creation and annihilation 
operators, a\ and a^, respectively. The Einstein oscillator 
frequency is w = (KjM) 1 ! 2 , and g = a/(2ifw ) 1 ' 2 - In 
addition we have introduced the electron kinetic energy, 
with nearest neighbour hopping amplitude t and chemi- 
cal potential fi. Finally, U describes the static electron- 
electron repulsion. 

Treating the four-fermion term involving the boson in 
mean field leads to a model with ordinary Holstein-likc 
electron-boson coupling 

Hei-b = g{n)uj {a\ + a l )(n lT + n a ) (3) 



9(n) = -g (4) 

where the coupling constant g(n) increases with electron 
occupation and is maximum for the Fermi level at the 
top of the band. Hence the quasiparticle dressing in this 
model increases as the Fermi level rises in the band, and 
conversely it decreases as the system is doped with holes. 
There are indications in the high T c cuprates that the 
quasiparticle dressing decreases with hole dopingjfj and 
this is one of the motivations to study this Hamiltonian. 
However to fully understand the physics of this model the 
mean field decoupling leading to Eqs. (3), (4) is certainly 
inappropriate. 

We study the Hamiltonian Eq. (2) in hole rather than 
electron representation. A particle-hole transformation 
yields in addition a Holstein-like coupling: 
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where now c. 



are hole creation operator and hole 



number operator respectively. One of the essential ingre- 
dients of this and other similar models is that the system 
is inherently not electron-hole symmetric. This is clear 
within the Hartree approximation, but the asymmetry is 
best revealed through a generalized Lang-Firsov trans- 
formation 0,0. 

In the following section we perform this transforma- 
tion, and derive a Migdal-like expansion valid in the anti- 
adiabatic limit, which allows us to study corrections due 
to finite boson frequency. While the antiadiabatic limit 
is very well understood [5j, the effect of non-infinite bo- 
son frequency (retardation) is not . The purpose of this 
paper is to address this question which is of fundamen- 
tal importance because derivation of the model from first 
principles necessarily leads to finite boson frequency 0- 
The expansion is performed to first order in inverse boson 
frequency, and the impact on superconductivity is deter- 
mined. We find that retardation increases T c at all hole 
densities, and increases the range of hole densities over 
which superconductivity occurs. This is in qualitative 
agreement with a previous study which determined 
the effective pairing interaction for two holes within a 
pseudospin model on small clusters. That study found 
an increased attractive interaction due to retardation as 
well. 



GENERALIZED LANG-FIRSOV 
TRANSFORMATION 

We use the transformation 

H = e G He- G , 
where Q 

G = ^g(a l - a|) (n^ + n H - n H n il)- ( 7 ) 
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Application of Eq. © results in 

H = uj ^2 a l a * ~ ^0^2 n%a + UcS E n *T n *l 
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where U c s = U — toog 2 , /xo = H + toog 2 , and 

Xl = cxp(g( ai -al)(l-n^S), CM 



dresses the hole-hopping amplitudes. These operators 
are dependent on both the oscillator and hole degrees of 
freedom. 

Eq. (JSJ) leads to a low energy Hamiltonian for the hole 
degrees of freedom if the ground state expectation value 



is taken with respect to the oscillator degrees of freedom. 
Following Ref. we find 
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Eq. UlUfl is the effective quasiparticle Hamiltonian stud- 
ied in Ref. |5j ■ This model is known to lead to supercon- 
ductivity at low hole dopings, and to characteristic prop- 
erties related to 'undressing' phenomenology 0. The 
effective single hole hopping 



t(n) = t o + nAt 
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is an increasing function of hole doping, and so is the 
effective hole bandwidth. 

Based on previous studies0, superconductivity (and 
in fact pair binding) occurs in the model Eq. (|SJ) in the 
anti-adiabatic limit. We wish to answer the question: 
will superconductivity occur more or less readily with re- 
tardation, i.e. away from the anti-adiabatic limit ? To 
this end, we use perturbation theory following the gener- 
alized Lang-Firsov transformation. That is we note that 
the operator xj can be written 



XL = ru 



+ (1 - rii-a) exp (g(a,i - a})) 
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and similarly for Xi a . In the last line of Eq. (|14|l we have 
expanded the exponential to first order in the oscillator 
momentum (or, equivalently, the displacement). To put 
terms involving the linear oscillator operators into 'stan- 
dard' form (with the oscillator displacement rather than 
the momentum) , we utilize the canonical transformation 
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Finally, Fourier transforming all the operators to momen- 
tum space, we obtain 



H 



9 



ho 



At 



to 



lAt 



N ^ 

kk' 



(16) 



3 



where 
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Note that we have included Hartree corrections already 
in the hopping and chemical potential terms; these are 
not to be included again when treating the Hamiltonian 
l|lfi|l . The correction to the hole hopping matrix element 
in particular leads to a wider bandwidth D for increasing 
number of holes in the band. 



It is also customary to write the normal self energy in 
terms of an odd and even (in Matsubara frequency) part: 

E(fc, iuj m ) = iw m (l - Z(k, iw m )) + x{k, ioj m ), (28) 

which leads to three coupled equations instead of the two 
Eqs. 123I24|1 . Note that the arguments of the g-factors 
in the kernels of Eqs. 123I24|I are actually different, and 
lead, in this case, to kernels with opposite signs in the 
pairing and normal equations. These equations also re- 
quire an auxiliary number equation: 



l + 2Rc 
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G (-k, -iw m ) 



k, m Go (k,iu m )G 1 (-k,-iuj m ) + (j) 2 (k,iuj m ) 



. (29) 



where n is the hole number density. 

These equations are most easily solved by noting that 
each unknown function can be decomposed into three k- 
dependent pieces. For example, 



ELIASHBERG-LIKE APPROXIMATION FOR 
HIGH BOSON FREQUENCIES 

Eq. nitty looks like a standard Hamiltonian with linear 
electron-phonon coupling. One can determine its prop- 
erties with the Migdal-Eliashberg approximation. Omit- 
ting the detailed steps in the derivation, the final result 
is: 
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where (f>(k, iu> m ) and E(fc, iu> m ) are the pairing and nor- 
mal self energies, and the denominator E(k, ioj m ) is given 
by 

E(k,iuj m ) = G (k,iu) m )G (-k, -w m ) + <j> 2 (k, iu m ) 

(25) 

with 

Go (k,iu) m ) = iw m - (4 - Mo) - E(k,iu m )- (26) 

These equations are standard Eliashberg equations, with 
iuj m = iirT(2m — 1) a fermion Matsubara frequency, /3 = 
l/(fcsT) the inverse temperature, and AT the number of 
lattice sites. The direct electron-electron interaction is 
given by 



Vkk' =U e _ 



2At 



t + nAt 



4 + 4' 



(27) 



(f)(k,iw m ) = 4>o(i^ m )+(f>i(iu) m ) 



(30) 

where D = D(n) — 8(t Q + nAt) is the n-dependent hole 
bandwidth. This makes the entire problem only slightly 
more difficult than the standard Eliashberg equations 
which require iterative solution in Matsubara frequency 
space 0. 

Nonetheless, we need not solve these rather compli- 
cated equations to determine their properties for large 
boson frequency, ujq. Eqs. n23!24|l both contain a kernel 
of the form: 



D(n) 



luo uji + vi. 2 



g 2 (to{l - n/2) + nAt/2 V 



+ nAt 



D(n)/2 



(31) 



where iv n = i2i:Tn is a Boson Matsubara frequency, and 
the other variables have been previously defined. We 
have written Eq. (f 3 1|) with an explicit prefactor D(n)/u)o, 
which shows that for infinite boson frequency, ojq, all of 
the complications due to the electron-boson coupling can 
be ignored. This leaves 

Z(k,iu> m ) = 1, 
X(k,iuj m ) = 0, 
cj){k,iuj m ) = (f>(k) 

i v- v <m 

— / *kk' i 1 • 

N P few G " (k', iuj m .)G ~ {-k', -iu m ,) + (t> 2 {k<) ' 

which is a straightforward BCS-like problem which has 
been solved previously The result is an extended s- 
wave solution for superconductivity at low hole densities. 
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To understand what happens for large but non-infinite 
boson frequencies, we adopt a standard approximation 
for Eliashberg theory, i.e. in the kernel we assume that 
the boson frequency is much larger than the other energy 
scales, so that it becomes independent of Matsubara fre- 
quency. This allows us to neglect the normal channels, 
and focus on the gap equation 



Substitution of Eq. |33 into Eq. |2U leads to 3 homoge- 
neous eqns, so that T c is given by setting the determinant 
of the following 3X3 matrix to zero: 

1 + «(n)/2 + uIq — kli n(n)l3 + ul\ — kl 2 k(ji)I^ + ul 2 — kl^ 
-kl - 2«(n)7 1 1 - kli - 2n{n)I 2 -kl 2 - 2n(n)h 



4>{k) 
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~Nf3 



D(n) 



where 



e/c' 



D{n)/2 J D(n) 



cj>{k') 



(33) 



n{n)I 

D(n)/2 
-D(n)/2 



dl 



D(n)/2 



n(n)Ii 



l-2/(e-/io) 



2 e 



Mo J 



1 + n(n)I 2 
(38) 



(39) 



where 



n(n) 



_ D(n) g 2 (t Q {l - n/2) +nAt/2V 



(34) 



uj 2 \ t a + nM 
The Matsubara sum in this equation is readily performed: 



E 



D(n)/2 _ D(n)/2 
+ (e/c - Mo) 2 2(i k - /t ) 



= Ji(e fc )- 



(1 - 2f(i k - Ao)) 
(35) 



As noted already in Eq. H30JI . the fc-dependence can be 
decomposed into a constant, a term proportional to e k , 
and a third term proportional to e 2 ,. The latter term is 
new, and arises from the fact that a non-infinite boson 
frequency is used. More explicitly, if we use a constant 
electron density of states (g(e) = l/D, —D/2 < e < 
D/2), then we have: 



<t>{x) 



dx' 4>(a/)Fi(?) 



(-K(n)(x 2 + x' 2 ) + 2n(n)xx' - k(x + x') - u cfi ^j ,(36) 

where u c g = U e g/D(n), k = 8 At/ D(n), and x = 
e/(D(n)/2), and similarly for the primed quantities. Nor- 
mally, positive components of the kernels contribute to 
pairing, and hence to T c . For example, focusing on the 
anti-adiabatic terms, w G g clearly deters pairing, whereas 
k will enhance pairing in the energy range where holes 
dominate (negative e and negative e'). At the same time 
k hurts pairing in the electron regime (positive e). Clearly 
the effect of n(n) is helpful (to pairing) in the central term 
(and hence adds to the role of k), and is detrimental in 
the first term. However, a full solution is required to de- 
termine the overall impact of a non-zero k(ti) (and hence 
non-infinite ujq — see Eq. (|34|) 'l flo| . 

BCS-LIKE SOLUTION 

A full solution to Eq. [JUlis obtained by noting that 

(f)(x) — ao — a\x + a 2 x 2 . (37) 



For the number equation we assume all non-pairing inter- 
actions are already taken into account in the parameters; 
thus we obtain, for the chemical potential: 



Mo 



D{n)/2 



= -(l-n), 



(40) 



where n is the hole density. The evaluation of the de- 
terminant, and the determination of T c {n) is a tedious 
process; in weak coupling, one can use the relations: 
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25 
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(41) 



where p = 1 — n. Weak coupling is accurate over a wide 
range of parameters . The end result is 

D{n) 



k B T c . = 1.13- 



■y/n(2 — n) exp (—a/b), 



(42) 



where 



a = 1 + kp(2 + kp) - Kp(-k - 2up + 3kp 2 ) + 

k 2 {-9 + SOp 2 + 7p 4 )/l2 - k 3 (3 - 9p 2 + 13p 4 + p 6 )/12 

b = (k 2 - 2u + Akp + k 2 p 2 )/2 + n{iu + 3up 2 - 4fcp 3 )/3 + 
k 2 (1 + 6p 2 + p A )/A - k 3 (9 - 9p 2 + 15p 4 + p 6 )/36, (43) 

with k = n(n). Note that if only the coupling to the bo- 
son were present, then, for small hole density, b/a « 2k 2 
for k « 1. Thus the boson leads to an attractive in- 
teraction, independent of the sign of k, but T c is unob- 
servably small if k << 1. However, if T c already ex- 
ists (through the presence of the anti-adiabatic terms, 
in this case), then the enhancement due to the coupling 
to the boson (which, in this case, means a coupling at 
non-infinite frequency) can be substantial. This nonlin- 
ear effect is merely due to the nature of the exponential 
function. 
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In Fig. 1 we plot T c vs. hole density n for parameter 
values At = 0.185 eV, U eS = 5.0 eV, t = 0.025 eV, with 
lo = oo,10 eV, and 5 eV. The explicit n dependence 
of /-c(n) is not significant, and is shown in Fig. 2, for a 
variety of parameter values (note: for Fig. 1, At /to = 
7.4) . Fig. 1 shows that T c is enhanced by retardation — 
as a->o decreases, kq = n(n = 0) increases. It is also clear 
that the range of hole densities increases with increasing 
amount of retardation. 

Finally, we can compute the effective mass and quasi- 
particle residue as a function of doping. These are ob- 
tained in a two-step process. First, in the anti-adiabatic 
limit (ujq —too), the effective mass ratio is given for low 
hole densities Ml as 



(1 + At/tp) 2 
(1 + nAt/to) ' 



(44) 



and the quasiparticle residue is similarly given by 



ZaaO 



(i + %At/t y 

(1 + At /t ) 2 



(45) 



Here z aa o is the quasiparticle residue on the Fermi sur- 
face, and the subscript aa stands for the anti-adiabatic 
limit. To incorporate changes due to retardation effects, 
we proceed in the usual fashion, and compute the real 
part of the analytic continuation to the real axis of the 
self energy expression Eq. I|24|) in the normal state, and 
without self-consistency (i.e. we omit the self energy on 
the right hand side). Then 



m * m le. (1-aw 



m m \ 1 + <2fc 



and 



where 



and 



zo 



^aaO 

I- a J 



(46) 



(47) 



a<j = -7^j-(e k ,u + i5)| w =o,6- h =#o (48) 



9Si 

an = - 5 ^(e~fe,w + i6)\u,=o,e k =iio) (49) 

where Ei is the real part of the self energy. These quan- 
tities are given by 



2/c(n) 1 + (1 - nf 



and 



K(n)[ 1 + (1 - nf 



(50) 



(51) 



The effective mass and residue are plotted in Fig. 3. It is 
clear that the effective mass decreases and the quasiparti- 
cle residue increases as a function of increasing hole con- 
centration, as expected in this model. Note that retarda- 
tion actually decreases the effective mass and increases 
the quasiparticle residue, because retardation 'undoes' 
some of the effects from the generalized Lang-Firsov 
transformation. Note that retardation has little effect on 
the doping dependence of these properties. To empha- 
size the observation made in Ref. 0] we have plotted the 
inverse effective mass ratio in Fig. 3b (dashed curve) in 
the anti-adiabatic limit, to show that the residue and in- 
verse mass behave similarly as a function of doping. This 
indicates that the frequency dependence of the hole self 
energy is most important (as opposed to the momentum 
dependence). 



CONCLUSIONS 

We have studied the effect of finite boson frequency in 
a dynamic Hubbard model. The model describes modu- 
lation of the Hubbard on-site repulsion by a boson degree 
of freedom. In the limit of infinite boson frequency the 
properties of the model are well understood, and here we 
studied the effect of finite frequencies through a gener- 
alized Lang-Firsov transformation and a high frequency 
expansion. 

One might also ask about the use of cruder approxima- 
tions. For example, one might first consider the Hamil- 
tonian (j^J in the mean field decoupling given by Eqs. 
(3,4). Doing then the usual Lang-Firsov transformation 
would yield an effective repulsive Hubbard model, with 
no superconductivity in the anti-adiabatic limit at the 
BCS-Eliashberg level of approximation. Thus the start- 
ing point for such a study would already be in serious 
error. The dynamics inherent in the coupling of the dou- 
ble occupancy of holes to the boson displacement given 
in Eq. JSJ is crucial for the occurrence of pairing. 

The central result of this work is the finding that T c 
increases with retardation . That this effect is expected 
was also indicated by the exact diagonalization study 
in Ref. [3|. In the derivation of this model from first 
principles^, the boson frequency represents the scale 
of intra-atomic electronic excitation energies. This scale 
is expected to be large compared to the scale of inter- 
atomic hopping, but is certainly not infinite. Therefore 
it is essential to study the effect of finite frequency cor- 
rections as done in this paper. 

This finding supports the possibility that the model 
may be relevant to the superconductivity of real mate- 
rials. It should also be noted that in the microscopic 
derivation of the model both the coupling constant g 
increases and the boson frequency decreases as the de- 
gree of negative charging of the ion increases. According 
to the results of this study, both effects enhance super- 
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conductivity in this model. Hence our results support 
the hypothesis that conduction of holes through negative 
ions is conducive to high temperature superconductivity. 
This is in qualitative agreement with the fact that high 
T c is found in cuprates (with holes in = ions) and in 
MgBi (with holes in B~ ions). 

We are grateful to the Natural Sciences and Engineer- 
ing Research Council (NSERC) of Canada and the Cana- 
dian Institute for Advanced Research (CIAR) for sup- 
port. 
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FIG. 1: T c vs. hole doping, n. For these parameters g = 
2.063, the hole bandwidth at the top of the band is D — 0.2 
eV and the electron bandwidth at the bottom of the band is 
D = 14 eV. The quasiparticle weight is 1 at the bottom of 
the band (electrons) and 0.014 at the top of the band (holes). 
The boson frequencies are uio = 10 and 5 eV, which gives 
«o = 0.043 and kq = 0.085 respectively. 
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FIG. 3: Effective mass (a) and quasiparticle residue (b) vs. 
hole doping, n, for several boson frequencies as indicated. 
As expected, the quasiparticle effective mass decreases as the 
situation becomes more electron-like; similarly the residue in- 
creases towards unity. Note that both quantities have the 
same tendency, as is made clear by the dashed curve in (b), 
where the inverse effective mass is plotted for loq — > oo. 



